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. Abstract 

^-j- ■ Simulating pp collisions at LHC energies in the framework of the SUGRA model 

and the detection of the produced leptons and jets by ATLAS we demonstrate that a 

Q\ ' clean signature of SUSY can be obtained over a large domain of the parameter space 

in the case of L-violating R parity breaking (A couplings). The obtained signal allows 
the reconstruction of the SUSY particles and thereby the precise determination of the 
Qh! model parameters too, mi/2, tan/3 and sign/^. 

1 Introduction 

^ ■ The ATLAS Collaboration has carried out a detailed study to detect the SUSY signature 

in the framework of one of the most popular model, SUGRA 0,0. It has been shown j|] 
that if SUSY exists at the electro-weak scale, it should be discovered by ATLAS and a 
general method has been given to determine approximately the mass scale of the SUSY 
particles. In subsequent papers it was shown in five representative points of the 

parameter space that some of the SUSY particles can be reconstructed and using the 
obtained characteristics (masses, branching ratios) the model parameters can be precisely 
determined || . All these studies have been carried out assuming that R parity is conserved. 

In this note we consider that R parity is broken in such a way that the lepton number 
L is violated through A- type couplings. The present experimental limits ||] cannot com- 
pletely exclude such a scenario. In this case one of the prominent signatures of SUSY, 
the missing energy is considerably weakened because the lightest SUSY particle (LSP) is 
allowed to decay. Due to this decay the lepton and/or jet multiplicity increases consid- 
erably and some efficient cuts (e.g. lepton veto against the t — t background) cannot be 
applied. On the other hand, the decay products of the LSP in some cases allow its direct 
reconstruction. Therefore the event topology and the search strategies are different of the 
case when R is conserved. This has motivated us to revisit the feasability to detect SUSY 
and to determine the parameters of the SUGRA model using the ATLAS detector. 
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In section 2 we give a brief description of the phenomenology of the R parity violation 
and the event generator used. Section 3 deals with the ATLAS detector and with the fast 
simulation of its response. In section 4 we present the domain of the parameter space 
where a SUSY signal can be expected by ATLAS. In the subsequent three sections the 
reconstruction of the SUSY particles and the determination of the model parameters are 
described in the LHC points 1,3 and 5 which represent a heavy, light and medium SUSY 
mass scale. We summarize the obtained results in the concluding section. 



2 Basic Phenomenology 
2.1 R parity violation 

i?-parity has been introduced |K| in order to avoid fast nucleon decay and flavor changing 
neutral currents (FCNC). If the multiplicative quantum number 

R = (_l)3B+L+25 (1) 

is conserved it guarantees automatically baryon number (B) and lepton number (L) con- 
servation. R is +1 for Standard Model (SM) particles and its value is -1 for their super- 
partners. The most important experimental consequences of the conservation of R are that 
super partners should be produced in pairs and the lightest superpartner (LSP) should be 
stable. The LSP interacts weakly, therefore the prominent signature of SUSY in case of 
R parity conservation is a considerable amount of missing (transverse) energy (E™ lss ). 

Although no violation of B or L has been observed yet, there is no firm theoretical 
argument which would require exact conservation of them and that of the R parity. In 
fact the following term in the superpotential 

Wfi = XijkLiL^ + KjkQiLM + KjkUtDpi (2) 

which violates explicitely B, L and R parity, cannot be ruled out experimentally. Here 
L and E are isodublet and isosinglet lepton, Q and D are isodublet and isosinglet quark 
superfields, the indices i, j and k run for the three lepton and quark families. The suffix 
c denotes charge conjugate. The first two terms violate explicitely L whilst the last one 
violates B. Present limits on the proton lifetime suggests that either the L or the B 
violating terms (i.e. the corresponding Ajjfc couplings) should vanish for the first family. 
Other experimental limits e.g. on lepton number violation: double (3 decay, ox on N — N 
oscillation, etc. indicate that the couplings in Equ. @ shouldn't be expected to exceed a 
few percent, and usually are much smaller than the gauge couplings. Even so, if R parity 
is violated the topology of the expected SUSY signal changes substantially. Since the LSP 
is no more stable, the missing energy is considerably reduced. On the other hand the 
decay products of the LSP increase the average number of jets and/or leptons in an event. 
In general, the event topology depends crucially on the size of the couplings. If, e.g. the 
couplings are of the order of ~ 1CU 2 or larger, the mass spectra, branching ratios, etc. will 
be different in the two cases where R is conserved or violated. If however the couplings 
are smaller than the above value, the dominant effect of the R parity violation is that the 
LSP becomes unstable. An estimation of the LSP lifetime as a function of the couplings 
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show [O that we can distinguish four subcases giving rise to different detection strategies: 

(i) 1(T 4 < A < 10- 2 

(ii) 1(T 6 < A < 10" 4 
(Hi) 10~ 9 < A < 1(T 6 
(iv) A < 10" 9 

In case (i) only the event topology changes w.r.t. the case of R conservation. In case 
(ii) one can observe a displaced vertex at the LHC energies. In case (Hi) the LSP decays 
outside of a typical LHC detector, however it can be catched by special purpose detec- 
tors [12]. Finally the case (iv) cannot be distinguished experimentally from the case of R 
parity conservation. 

In this study we have deliberately chosen to study case (i) and compare the result 
with the case of R parity conservation, because it represents a more difficult experimental 
situation than case (ii) where a displaced vertex could disentangle the LSP from the 
rest of the event. Moreover we have assumed that A- jfc = X'l- k = and only one of 
the Xijk coupling is different from zero in Equ. (||). Nonzero X' i j k and X'-j k are subject 
of other reports inside ATLAS [jl6|. The hierarchical structure observed for the Yukawa 
couplings in the SM motivates our hypothesis above. The Lagrangian corresponding to 
the superpotential of Equ. (||) can be written in terms of particle fields for our case as: 



L \ = 2 Xi i k eL i **Rk ~ eL * u % B *Rk + UL i ^ L 3 e Rk ~ e-U vl 3 e Rk + v u e Lj e Rk - e u v Lj e Rk ) + HC 

(3) 

where ul and bl,r are lepton fields, the tilde denotes the field of the superpartner, the c 
stands for charge conjugation, * for complex conjugation and i,j, k are the flavor indices. 
As stated before, if the A couplings are smaller than 10~ 2 , which is our case, the sparti- 
cle mass spectrum practically doesn't change and the main consequence of the R parity 
violation is the decay of the LSP. This process is depicted in Fig.[l| where we assume that 
the LSP is the lightest neutralino (x\)- The decay proceeds through an R conserving and 
an R violating vertex and in the final state there are always three leptons out of them at 
least two of different flavours, one neutral and the other two of opposite charges, since the 
LSP is supposed to be neutral. 

The prominent signature of this type of SUSY event is the spectacular increase of 
"stable" leptons: electrons and muons in the final state. The neutral leptons, neutrinos, 
give rise to some missing transverse energy, but its magnitude is much less than that if R 
parity is conserved. The flavour of the lepton in the final state depends on the values of the 
indices k. Since X^k is antisymmetric in i and j there are only 9 independent couplings 
which we choose as: A 121, A 122, A 123, A131, Ai 32 , Ai 33 , A 23 i, A232, A233- The first two families, 
1 and 2 give rise always to "stable" leptons, electrons and muons, (and the corresponding 
neutrinos) in the final state. If an index 3 appears, the lepton is not stable if it is a t, 
and its decay products are most of the time different from electrons or muons. Since e in 
Equ. @ is an isosinglet, e^ is a r. The number of the stable leptons is the most prominent 
for A121 and A 122- It is less spectacular, if an index 3 appears at the second place, and even 
less if the 3 appears in the third place. Finally, if two indices have value of 3 one has the 
least number of stable leptons. On the other hand, the number of the neutrinos, and with 
that the magnitude of the missing energy increases in the order of the above mentioned 
cases. The expected extra number of stable charged leptons and neutrinos per each LSP 
decay for the different A couplings calculated by the program [ 13 1 are given in Table g. It 
is clear that the average number of the leptons for different flavours gives a strong hint on 
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the coupling which is realized. E.g. the coupling A122 gives rise predominantly to muons 
whilst the coupling A123 results in equal number of electrons and muons, etc. 




Figure 1: The fl decays of y? (assumed as LSP) through A type couplings. 



Table 1: The number of produced stable leptons in the LSP decay if it is the Xi- We have assumed 
~ 18% branching ratio for the semileptonic decay of the r. 





Decay channel 


<N e > 


<N )1 > 


<JV V > 


121 


e ± u fl e T 




1.5 


0.5 


1 


122 






0.5 


1.5 


1 


123 




i/ e /i ± r =F 


0.68 


0.68 


2.36 


131 


e ± is T eT 


i/ e r ± e =F 


1.59 


0.09 


1.68 


132 


e ± iy T fi T 




0.59 


1.09 


1.68 


133 




i/ e r ± r =F 


0.77 


0.27 


3.04 


231 






1.09 


0.59 


1.68 


232 




^T ± /i =F 


0.09 


1.59 


1.68 


233 






0.27 


0.77 


3.04 



The presence or absence of r's influence how easily the SUSY particles can be recon- 
structed. On the one hand the signature is weakened by the taus since they produce less 
stable charged leptons, and these, partly originating from the r decay, do not produce a 
sharp endpoint for the Xi mass. On the other hand, taus produce more missing transverse 
energy, moreover the smaller number of stable charged leptons produced decreases the 
combinatorial background in the mass reconstruction. 
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2.2 The SUGRA Model 



The present study of R parity violation is carried out in the framework of the SUGRA 
model 0,0. Contrary to the minimal version of supersymetric models (MSSM) which 
has a very large number of unknown parameters, the SUGRA model is characterized only 
by 5 parameters which are the following: 

(1) mo, an universal scalar mass, 

(2) mi/2, an universal gaugino mass, 

(3) Ao a common trilinear interaction term, 

(4) tan(/3), the ratio of the vacuum expectation values of the two Higgs fields, 

(5) the sign of \x of the Higgsino mass parameter. 

The mass spectrum of the SYSY partners at the electro-weak scale as well as their decay 
branching ratios are obtained from the above parameters by solving the renormalization 



group equations (RGE). This is performed in our case by the program of IS A JET fl4 ]. 

The SUGRA model predicts a hierarchical structure of the masses of the SUSY par- 
ticles. The masses of the first two families of the squarks and of the sleptons are driven 
essentially by mo and m\/2 through an approximate relation : 

m h.R = m o + m } + c Ul,r) ■ m 2 1/2 + D{Jl,r) (4) 

where c(/) are some numerical factors of order 5.5 4- 6.0 for the squarks and of order 0.15 4- 
0.5 for the sleptons, D(/l r) are the so-called D-terms (in general less important). For the 
third family of the squarks and the sleptons there is a mixing due to the corresponding 
Yukawa couplings, pushing down the mass of one of the sfermions and the other one in the 
opposite direction. In some regions of the parameter space this could have the effect that 
the lighter stau (fi) could become the LSP but such scenarios are ruled out by cosmological 
considerations. The gauginos U(l), SU(2) and SU(3) are driven mainly by m^. Because 
the U(l) and SU(2) gauginos will mix with higgsinos to obtain the mass states (x^ , i = 1, 4, 
and Xj > J = 1; 2), part of that spectrum will depend on the \x parameter (in SUGRA fi is 
determined by the condition of electro- weak symmetry breaking and typically fi S> m^). 
We have the following approximate relations: 



m^o « 0.4 



2.4 • mi/a 

™>l/2 



m^o 
X3 



m^o ~ m-± « I u I 

a-4 X2 



(5) 



The Higgs sector is composed in SUGRA by five mass states (h°, A , H° and H^). In 
the case of the lightest one, m h o < 130 4 150 GeV in about all SUGRA parameter space. 
The masses of the other Higgses are in general very heavy and depend on tan(/3). 

The gaugino- higgsino mixing depends more strongly on the parameter values, and this 
in turn determines the decay branching ratios. 
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3 The analysis chain 



3.1 The event generators 

Simulation of the signal events 

Our basic program tool is ISAJET^] [l4j] which simulates for hadron colliders the pro- 
duction and decay of the supersymetric particles, as well as the underlying event, i.e. the 
accompanying partons and their hadronization. 

ISAJET does not include R parity violation in the SUGRA model. It has been intro- 
duced in detail in another event generator, SUSYGEN [18], written for e + e~ collisions. 
Therefore, inspired by SUSYGEN, a set of routines computing in detail all the fi (A-type) 



decays of gauginos (%' 



1,4 and xf,j 



1, 2) and sleptons was written for ISAJET 



This program is able to simulate production of supersymetric particles in a hadronic colli- 
sion with R parity broken, if this latter is accompanied by L number violation. To obtain 
correct predictions, the hypothesis of small values of the A couplings (A < 10~ 2 ) must be 
used, which is our present case, in order to neglect any correction in the sparticle mass 
spectrum brought by RGE's. The dominant effect of the R parity violation is the decay 
of the LSP. 

We have generated several million signal events. In Table § one can see their repartition 
in the SUGRA parameter space and for the type of the couplings. In the last column we 
quote the integrated luminosity in terms of LHC months the generated events correspond 
to where we have taken an LHC year equivalent to 10 4 /pb (low lumonosity run with appr. 
1/3 of efficiency). 

Table 2: The number of generated signal events 



m [GeV] 


m i/2 [GeV] 


A 


tan(/3) 


sga(n) 


coupling 


Events 


LHC month 


- 1500 


- 1500 





2 


+1 


Al23 


256 000 




- 500 


- 500 





2 


+1 


Al23 


300 000 




400 


400 





2 


+1 


Al22 


400 000 


~ 36 












Al23 


400 000 


~ 36 


200 


100 





2 


-1 


Al23 


1 000 000 


~ 1.5 












Al22 


1 000 000 


~ 1.5 


100 


300 


300 


2.1 


+1 


Al23 


450 000 


~ 36 












Al22 


450 000 


~ 36 



The samples in the first two raws were used to study inclusive reactions. The rest of 
the statistics is devoted to the 1st, 3rd and 5th of the so called LHC points, where the 
determination of the SUGRA parameters has been performed. 

Simulation of the background events 

As mentioned, the signature of the signal events is the appearence of a large number 
of high pt leptons accompanied by high pt jets and by a moderate amount of missing 
transverse energy. Such event topology is also produced, albeit in a much reduced level, 
by the decay of heavy SM particles and these events constitute the main background. 
We have studied such processes using the ISAJET and PYTHIA event generators 



1 We have used version 7.30 
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including initial and final state radiation. The background event statistics is listed in 
Table|. 

Table 3: The number of generated background events 



Reaction 


Number of events 


t-t 


1 200 000 


W - Z,W -W and Z - Z 


600 000 


Z-bb 


600 000 


Drell-Yan 


600 000 



3.2 Fast simulation of the ATLAS detector 

Due to the large number of events to be generated the detector response could not have 
been simulated in detail, using e.g. GEANT [pOfl . Instead a fast, so called particle level 
Monte Carlo program, ATLFAST ]15| has been usedQ In this program the ATLAS de- 
tector is described by a simplified geometry and apart of the acceptance the detector 
response is parametrized. Below we repeat the main features of these description. 

Description of ATLAS by ATLFAST 

The detector geometry is given in the variables of the pseudorapidity rj = —In tan(0/2) 
and azymuthal angle eft, where 9 is the solid angle of the particle produced in the interaction 
point. The granularity in 77 x <j) is 0.1 x 0.1 for \rj\ < 3 and 0.2 x 0.2 otherwise upto 
\rj\ = 5. The produced particles except the muons and the neutrinos deposite their energies, 
smeared by a resolution function, and are integrated in individual rj-cj) cells. The effect 
of the 2T solenoidal magnetic field on the deposited energy of the charged particles is 
parametrized. The effect of cracks in the calorimeter and in the trackers are taken into 
account by the parametrized acceptance function. 

Reconstruction algorithms 

First clusters are created from the cell energies using a simple algorithm. Next, isolated 
photons and electrons are reconstructed if the simulated particle falls in the acceptance 
region (typically \r/\ < 2.5), its energy, smeared by a parametrized resolution function, 
matches that of a cluster, and in the case of electrons, the r\ and (ft values of the cluster 
and the electron are the same within the resolution. The energy of the reconstructed 
particles is that of the simulated one smeared by a parametrized resolution, whereas their 
direction stays the same. 

After having removed the clusters of the isolated photons and electrons the recon- 
struction of the parton jets is carried out by a simple fixed cone algorithm. Isolation 
criteria of jets or charged tracks can be defined as a function of the deposited energy in 
a AR = \J (At/) 2 + (A</>) 2 cone. The reconstructed jet energy in general doesn't match 
with the true energy of the corresponding parton. Therefore a correction factor has to 
be applied in a later phase of the analysis. This correction factor, depending on the cone 
size and on the momentum of the jet, is established by comparing the reconstructed jet 
transverse momentum with that of the original parton: = pj parton /p^ et . It has been 

demonstrated that applying such factors one can correctly reconstruct the position of the 



mass peak of the Higgs boson 15]. This correction factor depends on the jet type and also 



We have used the version 1.57 
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if the jet contains prompt leptons. In Figj2| the correction is shown for the three different 
LHC points analysed. The observed difference in point 3 and the other two points is due 
to the facts that : 

(i) the production mechanisms of the b's are different. In point 1 and 5 they originate 
mainly from the decay hP — > bb, whereas in point 3 the production of the b's is less 
correlated; 

(ii) in point 3 there are more leptons produced and therefore the probability is higher 
that a lepton is inside the jet-cone. 

One has to stress that the dispersion of the correction factor can be rather large, reaching 
even 25%, especially at low pt (< 100 GeV). Fig^J shows the invariant mass of the bb pair 
from the h° mass at the LHC points 1 and 5 after having applied the correction. 





M(h°) = M'"(bb) in Point 1 M(h°) = M'~(bb) in Point 5 



Figure 3: The effect of the resolution and calibration of b jets on the invariant mass distributions 
in SUGRA points 1 (left) and 5 (right) and for the ft couplings X 122 = 10~ 3 : M(h°) ~ M(bb) - 
for reconstructed b jets originating from h 0, s. The ISR, FSR and hadronization are switched on. 
The arrows point to the theoretical values of m^o in each SUGRA point. 

Muons are reconstructed if they are in the acceptance region (typically < 2.5) 
and their energies and directions are obtained from the simulated one by a parametrized 
smearing function of ATLFAST. In most cases these functions were determined or cross 
checked by detailed simulation using GEANT. Finally, the missing transverse energy is 
calculated. 
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3.3 Detection and identification efficiencies 



The output of ATLFAST is written in a coloumn-wise ntuple (CWN) for further analysis 
using PAW |22]. This ntuple contains all the reconstructed electron, muon and jet objects 
together with the two transverse components of the missing energy. In addition, we have 
stored among the originally simulated particles and partons by RPV_ISAJET those which 
were produced in the decay chain of the generated SUSY particles. We have installed a 
bidirectional pointer between the reconstructed objects and the original particles/partons, 
and we have reinstalled the mother-daughter relationships between the stored original 
particles/partons. A special record containing the integrated luminosity has been also 
included. 

As a first step in the analysis we have randomly rejected reconstructed electrons and 
muons using a tabulated detection (including identification) efficiency. The efficiencies, 
as a function of the transverse momentum and 77, have been extracted from a version of 
ATLFAST (2.0) which was released after the bulk of our event simulation has been carried 
out (see Fig.£|). 



Electron and muon ID efficiencies vs. 77 and p, at 4> =0. 

£ id £ 
1 



0.8 - 




0.6 



0.4 - 



0.2 



p,*=60 GeV 



p,"=20 GeV 



p t '=40 GeV 



fc ID 

1 

0.99 
0.98 
0.97 
0.96 -_ 
0.95 
0.94 
0.93 
0.92 



: p t "=20 GeV 



VS. 7] 



-2 



p,"=60 GeV 




p,"=40 GeV 



£ id VS. Tj 




- ID 

1 

0.99 
0.98 -_ 
0.97 
0.96 
0.95 
0.94 
0.93 
0.92 



100 150 200 

P. (GeV) 



77=1.4 



£ ID VS. p t 



50 100 150 200 

P, (GeV) 



Figure 4: Detector efficiencies for muons and electrons as a function of pt and 77. The dependence 
on ip is practically negligible. 
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We have also randomly reclassified tagged b-jets to light jets with a b-tagging efficiency 
of 60% and vice versa using the misstaging probability of ~2%. These numbers have been 
obtained in a separate study by detailed Monte Carlo simulation 23]. 

The obtained final state reconstructed particles were submitted to selection cuts in 
order to find an optimum signal of SUSY over the background. These criteria as well as 
the reconstruction algorithms are described in the forthcoming sections. 



4 Inclusive measurements 



By measuring global variables, like e.g. the number of leptons of a given flavour or the 
average pt of a lepton in an event, or simply the number of events passing some selection 
criteria, we would like to answer the following three questions: 

(i) what is the maximum domain in the SUGRA parameter space in which ATLAS is 
sensitive for R parity violation; 

(ii) can we determine the approximate energy scale of a SUGRA signal; 

(Hi) can we determine the dominant type of the different couplings which causes the signal. 
For all these studies as well as for the exclusive measurements we have fixed the values of 
the Xijk couplings to 10~ 3 . We remind the reader that the event topology does not depend 
on the particular value of A if it is between 10 -2 and 10~ 4 . 

Sensitivity of ATLAS in the SUGRA space 

We have generated 1000 events in each of the 16x16 equally spaced points in the 
rrii/2 vs ttiq plane the other parameters being fixed (see Table^). We have generated the 
SM background events as given in Table|3[ Since the signature of the SUGRA event is 
characterized by multileptons and missing energy, we have tried several event selection 
based on the variation of the number of electrons and muons, their transverse energy and 
the value of the missing energy. We have used the quantity called significance: 

S = N stg /^I% (6) 

to delimite the sensitive region with S > 5, where N s i g is the number of the accepted 
SUGRA events and N^g is the number of accepted SM events. The sensitive region ex- 
tending to the highest value of mo and m\/2 has been obtained using the selection criteria: 

1. Ni > 3 

2. E r t niss > 100 GeV 

3. p\ > 70 GeV 

4. p[ > 20 GeV 

where Ni is the total number of e and /i, pi is the highest transverse momentum of the 
leptons, p\ is the transverse momentum of the other leptons in the event. The first two 
cuts reject almost all SM background coming from the decay of heavy SM particles. The 
Standard Model processes generally does not fulfill at the same time both constraints. Z 
pair production (by Drell-Yan or in the t, u channels) with Z — > l + + l~ give a higher 
multiplicity of leptons but with a very low EY 11 ^ (one of the leptons falling out the detec- 
tor acceptance could mime a but will decrease the lepton multiplicity). pair 
production (by Drell-Yan or in the t, u channels) and the associated production of W's 
with Z's give a higher i^™ 55 through the neutrinos but proportionally less leptons. The 
tt production with — > v\ + and the leptons arising from b jets reconstructed as 
isolated gives a higher E™ lss but is generally supressed by the isolation criteria of leptons. 
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The associated production Zjj is similar to the above mentioned cases. 



\ m = 10~ 3 , A = , tan/? = 2 , sgn/i = +1 




200 400 600 800 1000 1200 1400 



M (GeV) 

Figure 5: Sensitive domain in the SUGRA parameter space. For the explanation of the symbols 
see the text. The bricked domain is excluded by theory and the cross hatched one is excluded by 
experimental measurements. 

Fig.[| shows the domain of the sensitivity one can obtain in 1 year of running with 
LHC at low luminosity. Due to the limited number of simulated events (which is much 
less than one can detect in one year) we have an uncertainty on this region: full squares 
indicate the grid points where S > 5 at 99% CL, no symbols at the grid points indicate 
that S < 5 at 99% CL, finally the open squares correspond to the cases where one cannot 
make any of the two above statements. 

Energy scale of the SUGRA signal 

If a SUSY signal manifests itself, its energy scale can be determined in the case when 
R parity is conserved from the distribution of the quantity || : 

M ef f=i:pu+Er ss (7) 
i=i 

where pj. ; are the 4 jets with the highest transverse momentum. The missing transverse 
energy, £!™ ss originates mainly from the Xi decay. In the case of R parity violation we 
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should therefore replace E™ lss by the decay products of the x%i i- e - we have used instead 
of the definition (0) the following one: 

i=l i=l 

where p\ , are the 4 leptons with the highest transverse momentum. The distributions of 
M e ff for the LHC point No 1 and 5 can be seen in Fig|j. The point No 1 (3rd line in 
Table ^) is associated with a high mass scale, whereas the point No 5 (5th line in Table 
|2|) corresponds to a medium mass scale. This is well reflected in Fig|| 




Figure 6: Distribution of M e ff for the LHC points No 1 and 5 after 3 years of LHC run at low 
luminosity. The inclusive cuts used here for both points are (A122 = 10 -3 ): Ni eptons > 4 and 
j^miss -> 50G e y. The maximum in these distributions (M£fi x ) depends strongly on the mass 
parameters of the models being a good observable for the mass scale. 



In order to see how the M e ff distributions are correlated with the SUGRA mass 
scale, we have chosen randomly 30 points in the region of < mo < 500 GeV and 
< w-1/2 < 500 GeV. At each points we have generated 10 000 events (see 2nd line in 
Table |2|) and determined the maximum of the M e ff distributions: M^jJ x - The events 
were selected requiring more than 3 leptons and missing transverse energy higher than 50 
GeV in an event. In the first plot of Fig.^ we show the correlation of this quantity with 
the SUGRA mass scale, MsuGRA which we have defined as being the highest mass of the 
strongly interacting SUSY partners : 

MsuGRA = m.m(mg,mg R ,m- bi , m fTi ) (9) 
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The observed strong correlation is even more pronounced in the second plot of Fig.|7] where 
the distribution of the ratio of M^ff j MsuGRA is shown. 




Figure 7: Corellation between M^jJ x and Msugra and the distribution of the ratio 
M™ff /Msugra for 30 SUGRA points (for other explanations see the text). 



We have introduced another quantity to characterize the SUGRA energy scale in the 
case of R parity violation. Indeed, the algebraic sum of the lepton transverse momentum 
divided by the number of leptons: 

p l r rm = (tp\,)/N i (io) 

i=l 

gives a correlation with mo and mu2 as shown in Fig]8|, where the symbols delimite the 
regions where p l ^ norm \ s higher than a certain value. On this plot we have selected events 
with Ni ep t > 3, E^ 1133 > 50 GeV, and required a minimum value of 15 GeV for the 
momentum of any of the leptons and 50 GeV for the leading one. 

Sensitivity for the type of the coupling 

As mentioned in Section 2.1 counting the number of stable leptons not only can reveal 
the signal of R parity violation but also can give a hint on the type of the coupling which 
is realized. In order to show this we have simulated events in three different points of the 
SUGRA space (LHC No 1,3 and 5, see lines 3-5 of Table ||) and for two different couplings: 
A122 and A123. In the Table || the number of events are given one can observe for the two 
classes: Oe + 4// and 4e + 0/U in one year of LHC running at low luminosity. The events 
satisfy the following criteria: 

20 < p\ < 250 GeV and Ef lss > 250 GeV at points 1 and 5 

10 < v\ < 100 GeV and E™ iss > 150 GeV at point 3. 
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Figure 8: In the case of A123 = 10 3 , Aq = 0, tan/3 = 2, signa = +1 and for the global cuts : 
Ni ept > 3, E™ lss > 50 GeV, p\ > 15 GeV and p l t ' max > 50 GeV one can "divide" the (m ,m 1/2 ) 
space mp l t ' norm domains: - full squares: p t > norm g (30 ; 150) GeV; - fat squares: p l t ' norm g (150, 250) 
GeV; - open squares: j^' norm g (250, 500) GeV. Bricked domain is excluded by theory and the 
hatched one is excluded by the present experimental limits. 



Table 4: Number of events with different number of electrons and muons 





A122 


A123 


Oe + 4^ 


4e + 0/i 


Oe + Afj, 


4e + 0^ 


Point No 1 


934 ± 16 


3± 1 


76 ±4 


68 ±4 


Point No 5 


3010 ± 52 


20 ±4 


289 ± 17 


304 ± 18 


Point No 3 


52636 ± 1760 


3356 ± 444 


12400 ± 600 


15565 ± 672 



It is obvious that one can clearly distinguish which of the two couplings are realized in 
Nature. This is further illustrated in Fig|| where we have plotted the distribution of the 
number of muons/event for the two different couplings where the total number of electrons 
and muons is equal to four. 
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To conclude this section we have demonstrated that even if R parity is violated ATLAS 
can detect the signal of SUGRA in a large domain of the parameter space. This domain 
is compatible with that obtained in the case when R parity is conserved. We can equally 
well establish the mass scale of the SUGRA compared to the case of R parity conservation. 
Finally, one has a possibility to determine the type of the coupling if only one causes the 
violation of the R parity. 
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Figure 9: The distribution of the number of muons per event if the total number of electrons and 
muons is equal to 4 in the case of two /? couplings (A122, A123) in the SUGRA points 1, 3 and 5. 
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5 Exclusive measurements 



In this chapter we show that one can determine the parameter values of the SUGRA 
model, similarly to the case when R parity is conserved The quantities to be used 

for this purpose, in general, are the masses of the reconstructed SUSY particles as well as 
their observed production cross sections and branching ratios. We shall use only the first 
type of characteristics in this analysis since determination of cross sections and branching 
ratios is more sensitive to the acceptances and in many cases detailed simulations are 
needed. In 3 out of the 5 LHC points, which are sufficiently different to illustrate the 
methods in various conditions, we shall first show how to reconstruct the SUSY particles, 
determine their masses, and finally fit the model parameters to these mass values. 

The reconstruction of the SUSY particles are difficult because there are always at least 
two final state particles which cannot be detected. These are obviousely the LSP's in the 
case of R parity conservation, but if R parity is violated through the terms ([|) there is 
always at least one undetectable neutrino produced in each LSP decay (see also Table [[]). 
To overcome this difficulty we use the fact that in a 3-body decay: 

A — >abc (11) 

the invariant mass m^c of two of the three final state particles, e.g. b and c gives a clear 
endpoint. The endpoint, mf^ d is related to the masses of A and that of the undetected 
particle, a: 

mtf = m A -m a (12) 

since in the restframe of A the three-momentum of a (or that of the (be) system) is zero. 
This equation is particularly useful if the mass of a is zero: in this case one can directly 
estimate the mass of the particle A by measuring the endpoint of the distribution. Of 
course, in the practice this measurement is difficult because the endpoint may be hidden 
by eventual background and it can disappear if either b or c is unstable. This is the case 
when a r is produced in the decay of the x§, thus in all couplings when a "3" appears in 
the index of the A. 

Once the endpoint is established one can determine even the four-momentum P^ of 
the undetected particle a or that of the mother particle A by selecting events around the 
endpoint: 

Pia = —(K + Pn (M = l,-,4) (13) 
mbc 



As one can see from Equs. (^) and ([13]), if particle a has zero mass (as e.g. in the case 
of a neutrino), the four- momentum of the charged lepton pair is equal to that of particle 
A. The four-momentum of particle A can be used further for the reconstruction of the 
parent particle, and so on, until one arrives at the beginning of the decay chain, i.e. at 
the original SUSY particle. It is also obvious that at the endpoint the 4-momentum of 
particle a and A are parallel whatever their masses are. 



Since equations (12) and ( [13] ) are strictly valid only at the endpoint, one usually needs 
a large number of produced events and in the selection of the size of the region around the 
endpoint an optimum has to be found between the statistical error and the approximate 
validity of the above relations. In practice one applies in Eq.( |l3|) m(, c instead of m™ d in 
order to "rescale" the values of P^ + P^ which are smaller than their corresponding values 
at the endpoint. 
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If the decay (PH) proceeds through a sequence of 2 two-body decays: 

A — > Bb B — >ac 
one can observe an endpoint in the mt> c distribution, whose value is given by: 



(14) 



end 
m bc 



m n 



m B 



(15) 



The quantity of a x BR of the produced particle A helps to disentangle which of the two 



decay modes (11) or (|Tj) has occured and to choose between (12) or (15) to estimate the 
masses. 

At LHC the generic production and decay chains are represented in the schemas here 
below: (*). The squark or gluino undergoes the cascade decay: 



9 + Q 



± + q" 



U X? + W ± 



+ r + v 



Ul + + i- 



-)/h° 

+ v 



(*) 



and the R parity violation occurs at the end of the decay chain in the decay of the Xi- 

For the exclusive measurements we have chosen 3 out of the 5 LHC points to study. 
Point 1 has a high, point 3 has a low and point 5 has a medium mass scale. Another 
crucial difference between these points, coming from the mass parameters mo, m\i2 and 
sign(/i), is reflected in the decay of x®- I* 1 point 3, X2 decays mainly through a three 
body decay x% ~ * Xi + ' + ^ (with a virtual Z or I). In points 1 and 5 the main decay 
(65% -T- 80%) of X2 proceeds through a two-body decay X2 — * Xi + Du t also with a 
non-negligible branching ratio (20% 30%) it decays in a sequence of 2 two-body decays 
X2 -* Ir + lT -> Xi + ^ + l T with a being real. 

Therefore point 3 has a completely different event topology, namely a very large number 
of leptons in the final state, in comparison with the other two points. This is the reason 
that we present here an analysis in point 3 although it is already excluded by the LEP limit 
on the Higgs mass p4|] . In a nearby point, which is not yet excluded all our conclusions 
can be considered as valid. 

At each LHC points we have considered two different couplings: A122 and A123. In 
the case of the first one the reconstruction of the SUSY particles is easier than in the 
second one, where a r particle appears always among the decay products. This means 
more missing energy in the second case and also the absence of a clear endpoint in the 
invariant mass of the opposite-sign and different flavor (OSDF) lepton pair. This in turn 
makes very difficult if not impossible the reconstruction of the x? • 

The values of the SUGRA parameters as well as the number of generated events is listed 
in Tabld2|. The masses of the SUSY particles in the three points are listed in Table 
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5.1 The LHC points Nol and No5 



Concerning the reconstruction of the SUSY particles these two points are very similar and 
therefore we treat them here together. In the cascade decay chain (*) the predominant 
decay of the x% proceeds via x\ — * ^° + Xi with a BR of ~ 99% in point 1 and ~ 63% 
in point 5, where the competitive decay through In becomes also important. 

The main inclusive cuts to select SUSY events are based on the requirement of a high 
lepton multiplicity, Ni > 4 and a moderate missing transverse energy, ii"™ 55 . 

The SM background is small and can be completely neglected applying a cut on 
as is shown in Fig. 10 where the E^ 1133 distribution is shown for the points 1 and 5 and for 
A122 and Ai23- 



Table 5: SUSY particle masses at points 1, 3 and 5 in GeV 



LHC point 


1 


3 


5 


9 


1008 


299 


769 


QL 


958 


317 


687 


QR 


925 


312 


664 


62 


922 


313 


662 




855 


278 


634 




913 


325 


706 


h 


649 


260 


494 


II 


490 


216 


239 


Ir 


430 


207 


157 




486 


207 


230 


T2 


490 


216 


239 


h 


430 


206 


157 


xi 


775 


274 


526 


xt 


326 


96 


232 


xi 


778 


275 


529 


x» 3 


762 


258 


505 


X'2 


326 


97 


233 


Xif = LSP 


168 


45 


122 


h 


98 


69 


95 



The "structures" seen in the hatched histograms of Fig. 10 is due to statistical fluc- 
tuations. In what follows we shall demonstrate that in the case of the coupling A122 we 
can reconstruct all the SUSY particles in the decay chain (*), starting with the Xi> and 
in this respect we can achieve more than it was the case with R parity conserved. On the 
other hand, if A123 is nonzero, we can reconstruct only a fraction of the SUSY particles, 
and the decay products of the \ii i- e - the additional leptons, increase the background 
and thereby deteriorate the determination of some particle masses. We shall show that in 
spite of these difficulties the achieved precision is at worst comparable with the one in the 
case of conserved R parity. 
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Figure 10: The distribution of £™ ISS for SUSY events (open histogram) and SM events (hatched 
histogram) in points 1 and 5 for A122 = 10~ 3 (left plots) and A123 = (right plots) for events 

with Ni e ptons > 4, p l t ept > 10 GeV and after 3 years of LHC run at low luminosity. The arrows 
precise the cuts to be used. 



The case Ai 22 7^ 

Reconstruction of Xi u e(fi) + e± (/ u=t ) + M T 

For the reconstruction of the Xi we a Pply the following selection criteria: 

(i) Ni>4 (and = Nf and N e = for point 5), 

(it) p[ > 10 GeV, cos(a z ±,' T ) > 0.5, 

(in) E^ iss > 50 GeV , 
where « ; ± ; / T is the angle between any OSDF (i.e. electron -muon) leptons. 
The invariant mass distributions of the OSDF lepton pairs are shown in Fig. 11 for point 



1 and in Fig. 12 for point 5, where the number of events correspond to 3 years of LHC run 
at low luminosity. 

In point 1, isolated leptons are produced practically only in the two fi decays of the Xi s i 
therefore the background of OSDF distribution will be mainly Xi combinatorial. In point 
5, as we mentioned before, a non-negligible production of isolated leptons (beyond the 
ft decays of x?' s ) comes from X2 ~* + ^ T ~~ * Xi + ^ + This decay will produce 
two leptons of opposite sign (OS) but of the same flavor (SF), changing the balance of 
leptons per flavor. Most of the events with an increased number of leptons, will have only 
one X2 decaying leptonically (the other one, if it exists, will decay in h 0, s). To decrease 
the number of bad combinations (with leptons not coming from Xi' s )> additional cuts like 
= and N e = in each event are very efficient. 

There is a clear endpoint over a moderate background at both points 1 and 5 corre- 
sponding to the Xi mass (c.f. Table[|) in virtue of Equ.(12). We parametrize and subtract 



the background with a Maxwellian distribution and fit the resulting distribution near the 
endpoint obtaining the Xi mass values: 

m meas = 169 . 80 +0-2 GeV at point 1 (16) 

m ~ = 122.62i?'£ GeV at point 5 (17) 
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M(l*l") S4-B and B 



Maxwelllon fit of the background 




M(x°i) = 1 69. SO ± 0.07 GeV 



Figure 11: The invariant mass distribution of the OSDF lepton pairs for 3 years of LHC run 
at low luminosity at point 1 (for the selection criteria see text). In the upper- left plot are rep- 
resented all events (open histogram) and the background (shadded histogram) which is mainly 
SUSY combinatorial. The upper-right and lower-left plots show a Maxwellian-like distribution 
fitted to the background. The lower-right plot represents the result after the subtraction of the 
fitted background. The edge is fitted with a polinomial function and the error is only statistical. 



M(l*|-") S-l-B and B 



Maxwellion fit of the background 




50 1 OO 1 50 200 
M(ri ") S+B and B fit 




50 1 OO 1 50 200 250 
M(ri-~) after B subtraction 




1 OO 1 50 200 2 50 

122.62 ± 0.35 GeV 



M(x°.) 

Figure 12: The invariant mass distribution of the OSDF lepton pairs for 3 years of LHC run at 
low luminosity at point 5 (for the selection criteria see text). 
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The estimated error contains both statistical and systematical. This latter is due to 
the energy resolution of the leptons and mainly to the finite bin-size of our histograms 
around the endpoint. Since due to the bin-size one tends to overestimate the mass value 
the systematical error is asymetrical. 

For any further reconstruction using OSDF pairs (like X21 Xi > e tc) one W1 U take as Xi 
candidates the OSDF pairs with an invariant mass in (m^o — Am-^o, m^o), where Am^o = 
50 GeV for point 1 and Am^o = 30 GeV for point 5. To improve the statistics in point 5 
we will not use the additional cuts (see (£)). 

Reconstruction of h° — > b + b 



To the global cuts (JV, > 4, P\ > 10 GeV and E r t niss > 50 GeV) we will add some 
specific cuts on b jets. Since the h° decays to bb pairs with appr. 88% BR its reconstruction 
proceeds by the selection of these pairs, adding thus the following detection criteria to the 
above ones: 

(iv) p\ > 30 GeV for point 1, 40 GeV for point 5 respectively, and p\ < 300 GeV 
for both points 

(v) cos(a fe g) > 0.4 (point 1) or 0.3 (point 5) 

where a b i is the angle between the b and b. The obtained invariant mass distributions of 
the bb pairs are shown in Fig.|l4| and Fig. 15. There is a clear mass peak corresponding to 
the h° particle (c.f. Tabled). 



Sparticle sources of b jets in Points 1 and 5 
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Figure 13: The sparticle mothers of b jets (in Pythia codes) in the points 1 and 5 (see text). 



As we can see from the Fig. 13, the b jets are mainly produced through h — > bb but 
also through g — > b^) + b —> Xm) + b + b or g — > t\ + t —> x?( 2 ) + t + * ( or ~> Xi + b + t) 
and the subsequent decay of t. In the case of conserved R this SUSY background (as well 
as the SM tt) is efficiently rejected with a veto on leptons (coming from the top decay). 
However, in the case of fi through fj couplings we cannot apply anymore such a cut (and 
a veto on additional leptons doesn't increase significantly the S/B ratio). Another source 
of background (albeit small) is the production of Z's (i.e. 1 2 — * 1 1 + Z° or Xi — > Xi + Z Q , 
i = 2, 3, 4 or xf ~ > X± + with Z° — > b + b) resulting in a peak very close to the h° one. 
This will give rise to an asymmetric h° peak with a higher width as can be seen in both 
points. A similar effect arises also from the calibration error of b jets. 
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After havig parametrized and subtracted this background we have fitted the mass peak 
with a Gaussian curve and obtained the values: 

m meas = QJ Q g ± Lfi Qe y &t pQint j ^ 
m meas = Q4 ? ± Q e y ^ g ^g) 

The estimated systematic error is mainly due to the uncertainties in the fr-jet correction 
factor shown in Fig.^. 



M(bb) S + B and B M(bb) after B subtraction 




M(h°) = 97. 08 ± 0.62 GeV 



Figure 14: The invariant mass distribution of the bb pairs for 3 years of LHC run at low luminosity 
at point 1. On the left plot are represented all the events and the background (shadowed). The 
background is fitted with a Maxwellian function and subtracted. The result is drawn in the right 
plot. The arrow points to the theoretical value of m^o. The peak is fitted with a gaussian. The 
written error contains only the statistical part. 



M(bb) S + B and B 



M(bb) after B subtraction 
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M(h°) = 94.7 ± 0.64 GeV 



Figure 15: The invariant mass distribution of the bb pairs for 3 years of LHC run at low luminosity 
at point 5. 



Reconstruction of Xq. —> h° + Xi 

The reconstructed h and Xi allows the reconstruction of the mother particle, the x[>; 
since this latter decays to the formers by 99% BR at point 1 and 63% BR at point 5. We 
will require the x.i candidates (OSDF lepton pairs) to be "close" to the endpoint : 
(ui) m{OSDF) e (m en * - 50, m endp ) GeV - in point 1 and, 
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m{OSDF) G (m end P - 30, m end 'P) GeV - in point 5, 
and the h° candidates (66 pairs) to be around the mass peak of h° : 

(vii) m{bb) G (m peak - 15, m peak + 15) GeV - for both points 1 and 5. 

Looking at the relations between %2> Xi an d h masses, one can expect in point 1 a higher 
boost of Xi ~ h° pair, comparing to point 5, therefore we demand : 

(viii) cos(a hQ ^o) > 0.7 (point 1) or 0.5 (point 5). 

Finally one gets the invariant mass distributions of the ZioXi pairs as shown in the Fig.lC. 
The background in this case is combinatorial on the one hand, and comes from SUSY 
production of ti pairs with : ti — > 6 + W + + 6 + W~ — > 6 + 6 + l + + l~ + V\ + v\, on the 
other hand. However, the strong correlation required separately between the leptons and 
between the b jets will supress it considerably. 



X° 2 X°i + h° reconstruction in Points 1 and 5 
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Figure 16: The invariant mass distribution of the hoXi pairs for 3 years of LHC run at low 
luminosity in points 1 and 5. The tail is due to bad combinations with the wrong x\- The arrows 
point to the theoretical values of m^o- The masses are determined by a gaussian fit around the 
peak. The error is only the statistical. 



The peak around the x\ mass value (c.f. Table ||) allows to estimate this latter by 
fitting a Gaussian curve on the peak: 

m rneas _ ^ 2 ± g Q e y &t point 1 ^ 
m meas = 230 7 ± 3 9 Q e y at point 5 ( 2 1) 



Reconstruction of xt ~ * Xi + 

The same technique allows the reconstruction of the lightest chargino (xf)- As one 
can see in the scheme (*) one reconstructs first the W^. The W reconstruction is carried 
out by combining light quark jet pairs using the additional selection criteria (to (i -f- Hi) 
and (vi)): 
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(ix) p{ > 100 GeV and p{ < 600 GeV (point 1) or 350 GeV (point 5) 

(x) cos(ajj) > 0.9 (point 1) or 0.87 (point 5) 

where a ; - 7 - is the angle between the jj pair. The obtained invariant mass distributions 



of the jj pairs are shown in Fig. 17. There is a clear mass peak at the place of the W 
mass. Selecting the W candidates around the W mass peak (±15 GeV) and using the Xi 
candidates close to the endpoint (see (vi)) we plot the invariant mass distribution of the 
X\W pairs in Fig]l7], by requiring that the Xi and the W be close in phase space: 

(xi) cos^^o) > 0.85 in both points 
where a w ^o is the angle between the x^W pair. 
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X°i + W reconstruction in Points 1 and 5 
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Figure 17: The invariant mass distribution of the jj pairs for 3 years of LHC run at low luminosity 
at point 1 and 5 (upper plots). Selecting W candidates arround W mass peak (±15 GeV - between 
the bars) , we combine them with the Xi candidates and obtain the distributions represented in the 
lower plots for point 1 and 5 respectively. The arrows point to the theoretical values for m -± . 



One can see a clear peak at the mass value of the lightest chargino at point 5 (c.f. 
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Table[|). At point 1, where the statistics is much more limited due to the heaviness of 
the chargino, the corresponding peak is less clear. We have determined the mass of the 
chargino by fitting a Gaussian around the mass peaks which gives: 

m meas = g 2g>2 ± g 5 Ge y &t pQint j ^ 

Xi 

m meas _ 232 2 ± 45 Q e y &t pQint 5> ^ 



Reconstruction of q~R —> Xi + J e * 

A competitive cascade corresponding to this decay in (*) is : 



( 9 — ► QR + Q ) 

Ux°i+q (**) 
U /± + + v 



Due to the fact that the squarks are very heavy, the jet produced in the qn decay will 
carry an important fraction of energy. To avoid high combinatorial background one can 
ask only for one very energetic light jet in the event. 

Applying the following additional selection criteria (to (i -r- Hi) and (vi)) for the light 
quark jets: 

(xii) One light jet with pi > 750 GeV (point 1) or 400 GeV (point 5) whose 
invariant mass with any other light jet is outside the 

W ± or Z° masses (±15 GeV); 

(xiii) cos(a^o) > 0.0 in both points 
where a^o is the angle between the jxi pair. 
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Figure 18: The invariant mass distribution of the jxi pairs for 3 years of LHC run at low 
luminosity at point 1 and at point 5. The vertical arrows point to the nominal values of masses 
for different sparticles. 

The obtained invariant mass distributions of the jx\ pairs are shown in Fig.|l^. The 
qn can be produced directly in the hard process or in the decay of g. In points 1 and 5 
the jet associated with q~R is very soft comparing to the jet associated with x± an d often 
the soft jet can fall in the cone of the reconstructed hard jet. This is the reason of the 
presence of a shoulder at the g mass value in the distributions of Fig.^. Another source 
of background is the decay qi — * xf + Q with xf ~> Xi + — * £i + v l + ^ anci 
undetected. The left shoulder at the t± mass at point 5 corresponds to the case where the 
3 jets of the i\ — > xf + b — > x? + 3 + f + b or h — > X? + 1 ~ > Xi + 3 ' + f + b decay chains are 
erroneosely combined into a single energetic one. One can see mass peaks at the values of 
the right handed squarks (c.f. Tabld5|). The obtained mass values will have large errors 
due to the superposition of all these effects : 

m meas = g32 ± 2Q Q e y &t pQint 1 ^ 
m meas = m2 ± u Q e y &t pQint g ^ 



Reconstruction of q~i —> xi + 3 e ^ in point 5 

For the point 5 we can do more due to the higher statistics. As we already mentioned 
in the case of q~R, the qi can decay with a branching fraction of about 31 % in \\ and a 
light jet. This jet is also very energetic. In addition to the cuts (i -r- viii) one will select 
X2 candidates from the distribution : 

(xiv) m(xih°) e (m^o - 40, m^o + 40) GeV 
and only one hard light-flavoured jet with : 

(xv) p\ et > 100 GeV whose invariant mass with any other light jet is ±15 GeV 
outside of the or masses. 
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Figure 19: The invariant mass distribution of the jx% pairs for 3 years of LHC run at low 
luminosity at point 5. Vertical arrows point to the nominal values of qh and t\ masses. For cuts 
see text. 

The decay t\ — > x§ + 1 — ► X2 + 3 + J + & when the 3 jets are erroneusely combined into 
a single hard jet, manifests itself as a peak at the t\ mass (see Fig. 19). The gaussian fits 
to these peaks, after three years of LHC run at low luminosity, result in : 



m meas _ gg 5 ± 2Q Ge y &t point g 
m meas _ 5Q4 ± 2Q Ge y &t pQint 5> 



(26) 

(27) 



Reconstruction of 



Xi + I in point 5 



In the point 5, the l R decays ~ 100 % in Xi an d I ■ In order to obtain the mass of 
the right handed slepton, we have combined the Xi candidates with a "soft" lepton not 
participating in the Xi reconstruction and satisfying: 
(xvi) 10 < p[ < 200 GeV 



[xvii) cos(a^o y 



> 0.5 



where a^o is the angle between the lx\ pair. The obtained invariant mass distribution 



is shown in Fig.|20|. Besides the combinatorial one, the background arises from the decay 
Xi — » Xi + — ¥ Xi + ^ + V U which results in an endpoint in the l^Xi distribution at 
the mass of X\ ■ 



N„/Bin 



SIr -> X°> + I (100% BR) in Point 5 



si r . 



i 



Point 5 



<m(SI„)> = 156.79 ± 1.65 GeV 
u(Sl s ) = 15.59 GeV 
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Figure 20: The invariant mass distribution of the Ix® pairs for 3 years of LHC run at low luminosity 
at point 5. The arrows point to the nominal values of mj B and m~± 
A gaussian fit to the mass peak gives : 
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m rneas = 15 g_ g ± j g Q e y &t pQint 5 ( 2 g) 

which agrees with the value of the right handed slepton (c.f. Table|5|). 



The case A123 7^ 

As stated earlier, in this case there is always a r particle among the decay products 
of the Xi, and this spoils the endpoint in the OSDF lepton pair mass distribution which 
would permit us to reconstruct Xi an d an the other sparticles further (see Fig.|2l] and 
compare with the upper left plot of Fig.[T^). Therefore, the strategy in the case of ft, cou- 
plings implying the third family of leptons must be changed, and consequently we shall 
abandon the direct reconstruction of Xx> an d return to the strategy developed in the case of 
conserved R parity. In other words, we shall discard the decay products of Xi i n the first 
stage of the reconstruction. 



OSDF invariant mass distibution in Point 5 (X l23 =10 3 ) 



1200 - 




250 



M(OSDF) GeV 

Figure 21: The invariant mass distribution of OSDF leptons in the case of R coupling A123 = 10~ 3 
for 3 years of LHC run at low luminosity at point 5. Open histogram - all events; shadded histogram 
- background and cross hatched histogram -signal. The arrow points to the nominal value of x% 
mass. The presence of a r lepton in the \\ decay will spoil the endpoint at the mass of Xi- 



Reconstruction of X2 ~ > + ^ T — * Xi + ^ + ^ (Xi ~~ * v e{n) + M ± ( e± ) + r± ) chain 
The decay of the X2 produces in this case at least three isolated leptons. With these 
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leptons one can form pairs of opposite sign (OS) and different flavors (DF) - typical for 
Xi decay, or of the same flavor (SF) - typical for X2- As one mentioned before, X2 wm 
produce firstly an OSSF pair of leptons through a double 2-body decay and in a second 
stage, Xi wm produce again one or two leptons. In the conserved R parity case one has 
used the OSSF combinations. Typically, the lepton associated with Ir has a higher p t 
than the lepton associated with X? as h can be seen from Fig.^. In these conditions, 
it is useful to impose, beyond the usual cuts (i.e. lepton multiplicity, £ , ™ ss , etc), a 
minimal difference between the pt of the two leptons taken in the OSSF pairs. To decrese 
the number of bad combinations with leptons coming from the subsequent r decay, one 
increases slightly the p^ m cut on leptons. The OSSF leptons, coming from different 
decays, are less correlated therefore we relax the cut in angle between them. To eliminate 
the events with x\ ~ * X? + h° we demand furthermore that no h° has to be reconstructed 
in the event. 
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Figure 22: The pt distributions of the leptons associated with Ir or xi in the X2 decay for 3 years 
of LHC run at low luminosity at point 5. 



In Fig, 23, are presented the invariant mass distributions of OSSF lepton pairs, where 
the events were selected using the following critera: 

(») Ni > 4, 

(it) p\ > 15 GeV , cos(ai±^) > 0, 
(Hi) E™ iss > 50 GeV, 

(iv) Ap[ = Ipl 1 -p\ 2 \> 60 GeV, 

(v) No bb pair with invariant mass in (m^o — 15, m^o + 15) GeV. 

The huge background comes from the bad combinations with the the leptons produced 
in the x? decays. Since the X2 decays into the OSSF lepton pair via two 2-body decays, 
the observed endpoint in the OSSF distribution is determined by the relation of Eq.( [i~5| ) 
between m^o, mr and m. o. After the fit of the Maxwellian background and subtraction, 
the value of the endpoint is estimated as : 



meas 
m OSSF 



111.9 ±2.5 GeV 



(29) 



where the error of about 2.5 GeV is dominated by the histogram bining (i.e. statistics). 
One can isolate lepton candidates coming from the double 2-body decay of X2 by a 
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selection around this endpoint. With the remaining leptons one can attempt to form OS 
lepton pairs and further combining them with the selected OSSF pairs to reconstruct the 
X2- Because of the neutrinos (always present) the X2 reconstruction is not complete and 
therefore an endpoint will appear in this OSSF+OS lepton invariant mass distribution 
depending on the m^o. 



M(OSSF) S + B and B 



Maxwellian fit of the background 



-2.673 ± 
-.2687E+05± 345.0 
5430, ± 59-03 
.1042 ± .9290E-O2 
24.1 1 ± .9504 
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Figure 23: The invariant mass distribution of the OSSF lepton pairs for 3 years of LHC run at 
low luminosity at point 5. 



If one applies the additional cuts to (i v) : 

(vi) mossF £ ( m o5SF ~~ 25\ m 05SF) GeV and rescaled as explained in the intro- 
ductory part of subsection 5, 

(yii) OS pairs constructed with leptons having p\ > 10 GeV, Ap\ > 20 GeV and 
cos(aos) > 0, 

(viii) cos(a ssF,os) > 0.85 , 



one obtains the distribution represented in Fig.|24 
One can observe an endpoint at: 



(30) 



m,o = 228.2 ± 5 GeV 

by linear extrapolation of the distribution as is the case for a 5-body final state (see p5]| ). 
This endpoint corresponds to the kinematical limit of the decay, namely to in^o. The 
error of the extrapolation is reflected in the large value of the error in m^o . 

A.2 
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Reconstruction of 



X2 + Q 



i end P (Eq(fL5|)) which means that the three- 
momentum of the Xi is nearly zero in the X2 restframe. Therefore selecting OSSF lepton 



We observe that at point 5 m^o 



pairs near the enelpoint of the distribution of Fig. 23 one can determine P~ using Eq.(|i~3|) 

_ — *2 

with the value of m^o from Eq( [3C]). These 4- momenta are then combined with a hard 



X2 

light jet not coming from a reconstructed W or Z°. To the criteria (i 
following : 

(vi — b) mossF £ { m> ossF ~ 30> m os5ir) GeV and rescaled, 

(ix) p 3 t 

around W, 

(x) cos(a j0 ssF) > 0. . 



(v) one adds the 



? et > 400 GeV and no invariant mass with any other light jet in ± 15 GeV 



The resulting distribution is depicted in Fig. 25. By a gaussian fit on the peak one obtains 
the value: 



meas 

N._/Bin x°2 (^sr„+r) x°i + ossf -> v + os + ossf 

140 
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(31) 



Point 5 

(A, 23 = 10- 3 ) 
<m„ a > = 228.2±5 GeV 
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Figure 24: The invariant mass distribution of the OSSF+OS lepton pairs for 3 years of LHC 
run at low luminosity at point 5. Due to the presence of neutrinos (from \i decay) an endpoint 
depending on m^o and m^o appears (see text). 
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N„/Bin s 1i ^ X 2 + j» 
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Figure 25: The invariant mass distribution of the (jet OSSF) pairs for 3 years of LHC run at 
low luminosity at point 5. The arrow points to the nominal value of mass. The mean value is 
obtained by a local gaussian fit on the main peak. 



rriq L can also be determined by combining the (OSSF+OS) pairs selected near the 

whose 4-momenta are approximately equal to P_ according to 

X2 



m^o endpoint of Fig.| 



Eq.( 13), with a hard light jet not coming from a reconstructed W or Z°. For that one 
will use the cuts (i -r- ix) supplemented by: 

(2*) m^o G (m e ossF,OS ~ 50 > m OSSF,Os) GeV > 
(xii) cos(a j0 ssF~os) > . 
The invariant mass distribution is shown in Fig. 26. One obtains the value: 



meas 
" L qL 



686 ± 12 GeV. 



(32) 



The agreement between the results of ( |3lJ) and ( |32j) justifies the procedure that we have 
applied in obtaining the mass of 
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N _ /R , n Sq L -> X °2 + Wd - (X°2 = OSSF + OS) 
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Figure 26: The invariant mass distribution of the (jet OSSF-OS) pairs for 3 years of LHC run at 
low luminosity at point 5. The arrow points to the nominal value of ql mass. 



Reconstruction of h° — > b + b 

The reconstruction of the h can be carried out almost with the same precision as 
it was the case of the coupling X±22- However, the applied isolation criteria reject more 
events due to the finite size of the r jets as compared to the case of X±22- This can be 
seen in the bb invariant mass distributions which should be independent of the fi coupling 
type. Following the same procedure as in the case Ai22> with the selection criteria slightly 
modified ; 

(xiii) 15 GeV < p b t < 300 GeV, 

(xiv) cos(«6b) > 0.4 GeV, 

one obtains the distribution shown in Fig.^7]. After the subtraction of the Maxwellian 
background and the gaussian fit of the peak one gets: 

m™ eas = 94.3 ± 1.5 GeV (33) 

This number includes the systematic error, mainly due to the incertitude on the energy 
scale of 6 jets. 



Point 5 

<m(Sq L )> = 686.1 ±10.3 GeV 
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N. v /Bin N e./Bin 

" I 60 I" Constant 47.89 ± .7835 

Mean 94.26 ± .1681 
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Figure 27: The invariant mass distribution of the (bb) pairs for 3 years of LHC run at low 
luminosity at point 5. 



5.2 The LHC point No3 



The salient feature of this point w.r.t. the other two is that the SUGRA parameters mo 
and mi/2 are much lower, hence all sparticle masses are lighter and the production cross 
sections are highly increased. As a consequence we could generate only a fraction of the 
events one would be able to detect with the ATLAS detector in one year even with low 
luminosity. In addition, due to the sigrifj, = — 1 and m^o — m^o ~ 52 GeV, the second 
lightest neutralino decays into OS and SF lepton pairs and not into a Higgs (h°, with a mass 
~ 70 GeV). Therefore in this point one cannot determine the Higgs mass. The number 
of leptons in the final state being higher, the combinatorial background is considerably 
larger in the channels where the particle reconstruction involves leptons. Furthermore, as 
we precised before, Xi can decay (through fi A-type couplings) in lepton pairs of different 
flavors as well as of same flavors. The consequence is a much more complicated structure 
in the invariant mass distribution of the OSSF lepton pairs than in the case of conserved 
R parity. 

Another feature of point 3 is that in the decay chain (*) the dominant decay products 
of the gluino are the b\ and the b. Since in the decay of the b\ the gaugino (x° or Xi) 
is accompanied again by a 6-quark, there is a large number of 6-quarks produced in each 
event. 

In the case of conserved R parity @ the observables which can be used for the determi- 
nation of the SUGRA parameters are the functions m g (m^o ) , (m^o ) and the difference 
Am^.o = m^o — m^o (from the invariant mass distribution of the OS and SF lepton pairs: 
x\ — * (^ + ^~) + X? - see the decay chain (*) and Equ. fll2])). In the ft case, one can develop 
two strategies depending on the A type coupling : 
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(1) For A couplings not producing r jets in the LSP decays (i.e. A122) one will directly 
reconstruct the x? as i n the other LHC points. The consequence of this new information 
is a better fit of the SUGRA parameters comparing with the R conserved case; 

(2) For A couplings producing r jets in the LSP decays (i.e. A123) one returns to the 
strategy developed in the R conserved case. 

As in the LHC points 1 and 5, we will present distinctively the two representative 
cases: A122 and A123. 



The case A122 7^ 

The analysis is focused firstly on the leptons. With these, one can form OS pairs of 
the same flavor (SF) or different flavors (DF). Due to the facts that: 

- Xi has a significantly higher production rate than the x°> an< ^ 

- only X? produces OSDF lepton pairs, 

the OSDF pairs will tag much better their origine than the OSSF ones. Therefore, in the 
first stage one will reconstruct the invariant mass distribution of OSDF leptons. This gives 
the X? mass through the same type of endpoint structure as in the other LHC points. The 
main background in this distribution is combinatorial. One will select OSDF lepton pairs 
in the neighbourhood of this endpoint. With the remaining leptons (originating with a 
higher probability from X21 but also fom x?) one will form OS lepton pairs. The endpoint 
in their invariant mass distribution determines Am^o = m^o — m^o . Selecting OS pairs 
near this latter endpoint and combining them with the x? candidates (from the OSDF 
distribution) one can completely reconstruct X2- Afterwards, combining \\ candidates 
with one and two b jets one will reconstruct b\ and g. 

For this analysis one uses the following global cuts : 
(i) Ni > 4, and vT > 10 GeV, 
(it) E™ iss > 50 GeV. 



Reconstruction of x? — * v eUi) + / i± ( e± ) + A t± 

Using the global cuts (i ii) and : 
(in) cos(aosDF) > 0.85, 



one obtains the invariant mass distribution of OSDF leptons as in Fig.|28|, where the back- 
ground is mainly combinatorial. After the fit of the Maxwellian background, subtraction 
and the polinomial fit of the resulting endpoint one obtains the x? mass : 

m ~ = 44.8+0.1 GeV ( 34 ) 
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M(OSDF) S + B and B 



Maxwelllan fit of the background 



-6.391 ± .1632E-01 

.1008E+05 ± 6.749 

890.2 ± .6514 

.1047E-01 ± .1192E-03 

186.2 ± 2.473 
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Figure 28: The invariant mass distribution of the OSDF lepton pairs for 1 year of LHC run at 
low luminosity at point 3 in the case of fi coupling A122 = 10 -3 . 



Reconstruction of x\ ~ * Xi + ^ ± (e ± ) + /i =F (e =F ) 

The x? candidates are identified by requiring OSDF lepton pairs with : 

(iv) miOSDF) e (mT - 10,m™ eas ) GeV, 

with rescaled 4- momenta according to Eq,(|l3|). One then selects leptons produced directly 
from the 3-body decay of X2 am ong the remaining OS lepton pairs demanding that the 
angle between the leptons (aos) satisfies : 

(v) cos(aos) > 0.5 . 

The invariant mass distribution of these OS pairs (Fig.p9|) has a complex structure due 
to the fact that, beside the combinatorial background (Maxwellian type with a long tail) 
there is also a contribution from the Xi decay itself (in fact the contribution is 50% of its 
branching ratio and it has the shape comparable with that of OSDF distribution). There- 
fore the OS distribution shows a sharp edge at about 52 GeV, correspondig to the mass 
difference m^o — m^o and a second one, less pronounced around 45 GeV, corresponding to 
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the Xi decay products. The tail beyond 52 GeV is due to the combinatorial background 
and we have removed it. The polinomial fit of the higher endpoint gives the value: 

m ~ = m -o - m^o = 53.3 ± 0.6 GeV (35) 
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Figure 29: The invariant mass distribution of the remaining OS lepton pairs, after the removal 
of \i candidates. The distributions corresponds to 1 year of LHC run at low luminosity at point 
3 in the case of ft coupling A122 = 10~ 3 . For explanations see text. 

As a next step we apply the cuts (i 4- v) and : 

(ui) m{OS) € (m$f as - 10,m^| as ) GeV, 
in order to reconstruct the \% by combining the OSDF and the remainig OS lepton pairs. 
We require that the angle between the two pairs should be small : 

(vii) cos(a SDF,os) > 0.8. 



The obtained mass distribution is shown in the Fig. 30. The gaussian fit around the peak 
results in the value : 

(36) 



meas = gg ? ± Q 2 q y 
X2 



Out of the three quantities defined by the Eq.(p3), pq) and (pq) we use two independent 
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ones (Eq.(p3) and (BH)) , for the determination of the SUGRA parameters. 
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Figure 30: The invariant mass distribution OSDF + OS lcpton pairs. The distributions correspond 
to 1 year of LHC run at low luminosity at point 3. For explanations see text. 



Reconstruction of the chain g —>■ b\ + b —>■ X2 + b + b 

The configuration of the sparticle mass spectrum (see Table ||) and the decay branching 
ratios for that point, permit us to make the following remarks on the production of b jets: 

(1) practically all b jets originate from the decay chain of g (with a small fraction 
coming also from t\ and/or t), 

(2) the b jets from the first decay of g are very soft comparing to those produced in the 
second decay (and associated with X2), as it can be seen in the Fig, 31. Correspondingly 
we label the b jets as follows : 

{viii) hard b jets if : p\ > 50 GeV ; 
\ix) soft b jets if : 10 GeV < p\ < 50 GeV . 
Selecting X2 candidates around the X2 peak from Fig.|30| : 

(x) m{OSDF, OS) € (mP eak - 15, mP eak + 15) GeV, 

one performs a first reconstruction of the g mass as follows. First of all, one selects pairs 
of b jets with all b jets passing the cut pt > 10 GeV and in each pair we identify the hard 
and the soft jet. Next we combine the x[>> ^hard an d b so f t 4-momenta if : 

(xi) cosfa^oi ) > , and 

v A2 hard / 



(xii) cos^/^oi \ h . 



> 0.5. 



These cuts are justified by the distributions shown in Fig.|32p and |32| b 
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N b j et »/Bin P, of b jets in SUGRA Point 3 
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P t b (GeV) 

Figure 31: The pt distribution of b jets in point 3 for 1 year of LHC run at low luminosity after 
applying the global cuts N t > 4, E' t niss > 50 GeV and p\ > 10 GeV. One can clearely distinguish 
the b jets from g and from b±. 
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Figure 32: The first reconstruction of the gluino mass in point 3, for ft coupling A122 = 10 -3 and 
after 1 year of LHC run at low luminosity (plot c). Also shown are the angle distributions between 
the X2 and the hard 6-jets (plot a) and between the (x°, bhard) pairs and soft 6-jets (plot b). For 
further comments see text. 
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The resulting mass distribution is represented in Fig.^c. For the b\ reconstruction, 
one selects the x\ candidates according to (x) and the hard b jets (see (viii)) around the 
g peak of the FigJ3z 



\xiii) m; 



£ (m peak — 10, m peak + 10) GeV, requiring this time 



(xiv) cos(a^o hu ) > 0.5. 

x ' v A2 hard' 

The result is represented in Fig. 33. After the gaussian fit around the peak one obtains the 
value : 



in 



bi 



276.6 ± 3.0 GeV 



(37) 
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Figure 33: The invariant mass distribution of X2 an d bhard jets. The distribution corresponds to 
1 year of LHC run at low luminosity. For the selection criteria see text. 

This permits us to refine the reconstruction of the g by tagging it with the b\. We 
combine the X2 candidates (cut (xi)) with a hard b (cut (viii)) requiring an invariant mass 
within the 61 mass window: 

(xv) m^o b e (nC eas - 10, m^ eas + 10) GeV. 

v ' X2 u hard 01 01 

We add subsequently a soft b jet with the same angular correlation as (xii) to obtain the 



final gluino reconstruction shown in the Fig. 34. The gaussian fit around the peak gives 
the measured value : 

m rneas = j ± 3 q Q eV (33) 



40 




Figure 34: The final reconstruction of the g after having selected the (x®, bhard) pairs around the 
reconstructed b\ mass peak. The distribution corresponds to 1 year of LHC run at low luminosity. 



The case A123 7^ 

Reconstruction of \2 ~~ * X? + ^i^) + /• t=F ( e=F ) 

Comparing to the case A122 the presence of r jets in the decay products of Xi spoils 
the endpoint structure in the OSDF lepton pair distribution and therefore the direct 
reconstruction of the Xi- Although the electrons and muons from the r jets increase the 
combinatorial background of the OSSF lepton pairs, due the fact that they are softer than 
those from the 3-body decay of the X2 they do not spoil the endpoint structure. We select 
events with the usual cuts (i -r- ii) and OSSF pairs by requiring : 

(xvi) cos(oeosSF) > 0.5. 
After a fit of the Maxwellian tail, subtraction and the fit of the result with a polinomial 



(see Fig.^5J) , one obtains using Eq.fll^ 



™OSSF = ™% - m tl = 52.9l°i GeV . (39) 



Reconstruction of the chain g — > 61 + b ^ xi + b + b 

Since the A couplings do not affect the b jets (if we neglect the difference in the 
reconstruction and tagging efficiency due to the larger size of a r jet comparing to a 
lepton) the classification of b jets in soft and hard types is done by the same criteria as 
in the case Ai22- The reconstruction of this chain follows the same treatement as in the 
case A122 with the notable difference that this time one cannot directly reconstruct the X2- 
Instead, we have to assume a mass value for the x% m order to obtain its 4-momentum 
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at the endpoint of the OSSF lepton pair mass distribution according to the Eq.(|T^). We 



have taken m^o = 97 GeV. 

A2 
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Figure 35: The invariant mass distribution of the OSSF lepton pairs for 1 year of LHC run at low 
luminosity at point 3 in the case of /? coupling A123 = 1CP 3 . 

The cuts used in this case are the following : 

(x) * m(OSSF) € ( m> ossF ~ ^> m OSSi= 1 ) GeV, for the X2 candidates, 
and for the first reconstruction of the gluino : 

(xi) * co S (a^o bha J > 0.5, 

(xii) * cos(a^ot ) > 0.5. 



The invariant mass distribution of xfyhardbsoft pairs is shown in Fig.3£. Selecting the 
hard b jets if around the gluino peak : 

(xiii) * m- g G [mP eak - 10,m peak + 10) GeV and, 

(xiv) * coB(a^ 6hopd ) > 0.5. 

one obtains the 61 mass peak as shown in the Fig, 37. The gaussian fit around the peak 
gives : 

m^ eas {97) = 277. 5 ± 3.0 GeV (40) 
The 97 inside the brackets indicates that this value was obtained assuming m-o = 97 GeV. 
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Figure 36: The first reconstruction of the gluino in point 3, for ft coupling A123 = 10 3 and after 
1 year of LHC run at low luminosity . 
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Figure 37: The reconstruction of the b\ selecting the b jet around the g peak. The distribution 
corresponds to 1 year of LHC run at low luminosity. 
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Finally, selecting the hard b jets and the OSSF pairs around the b\ peak : 
(xv)* m^ b € ( m ~ - 10, mf eas + 10) GeV, 
we combine them with the soft b jets with the same angle cuts as above {(xi)* and xii)*). 



The resulting mass distribution is shown in the Fig. 38. The gaussian fit around the peak 
gives the measured value : 



mf ttS (97) = 301.1 ±3.5 GeV . 



(41) 



Sgl -> Sb, + b, oft -» x 2+b ted +b, o(t -» +b tart +b, 
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Figure 38: The final reconstruction of the g after having selected the OSSF lepton pairs and hard 
b jets around the 61 peak. The distribution corresponds to 1 year of LHC run at low luminosity. 

We have repeated the above analysis for several different values of rn^o- The result is 
shown in Fig. [353. Using a linear fit one obtains the following expressions : 



m 



"" "* f m**) = m^ eas (97) + %(m^o - 97) GeV (43) 



with 9~ bx = 1.48 ± 0.52 and 6~ g = 1.55 ± 0.55, repectively. 
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Figure 39: The dependence of mg i (full squares) and rrig (full circles) as functions of m^o. The 



arrow points to the nominal value of m^o 
by a linear fit. 



97 GeV. The factors 0j and 9g (see text) are extracted 



5.3 Determination of the SUGRA parameters 



In subsection 2.2 we have already reviewed the dependence of the sparticle masses as 
functions of the model parameters mo, mi/2, Aq, tan/3 and sign/i. Practically all sparticle 
masses are sensitive to m^- mo drives mainly the sfermion masses, sign// affects the 
gaugino branching ratios (i.e. X2 ~~ * X? + ^ + + l~ versus X2 ~^ Xi + ^°)> but also the 
sparticle mass spectrum. On the other hand, the observables are practically independent 
of Aq (c.f. Fig 42 and Table H). An obvious way to determine the above parameters 



consists of finding the minimum of : 

X 2 (m ,m 1/2 ,tan(3,A ) = 



ES=i[O^(m , m 1/2 , tan(3, Aq) - O^ as ]a^[Of(m o ,m 1/2 ,tan0,A o ) - 0% 



meas] 1 \r)thi 



(44) 



where (mo, m\/2-> tan/3, Aq) represents the theoretical evolution of the observable a in 
function of the SUGRA parameters, 0™ eas is the measured value of the same observable 
and o a p is the covariance matrix . 

3 In general 0™ eas are not independent, however at the present time we neglect the offdiagonal elements 
of a a 0- 
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The minimum of the \ 2 can be found by scanning through the entire parameter space 
II. The overall minimum is practically insensitive to Aq. The sign/i parameter can be 
determined unambiguosly ||. Once the overall minimum is found one can determine the 
error on the parameters at the minimum by taking the square root of the diagonal elements 
of the matrix A y - : 



Y? ■ -1 ■ i- J . oo!! 1 

a,P=l 



Ai i = [ 2^ Aafoap) , ^ai = -^pj, Pi = m, , m 1/2 , tcm/3, Aq, 6 § , 6~ h (45) 



Table |l^ summarizes all the measured values and errors of the chosen observables for 
all studied cases. For point 3 we have added the h° mass determined by an independent 
measurement (h° — > 7 7) with an estimated error of 1 GeV. As one can see from this Table 
in most of the cases the models are overconstrained : one disposes of more measured 
quantities than the number of parameters to be determined. Comparing the measured 
values of the observables with the theoretical ones (see Table ||) one can see that the 
errors are slightly overestimated. 

O th (pi) are the theoretical dependencies of the quantitities O meas collected in Table |lC| . 
However, there is an exception at point 3 (for A123 7^ 0) where, according to Eqs.(H5) and 



(42) one has : 



Of (pi) = m~ g (pi) - 0§(m%o(Pi) - 97), pi = m ,m 1/2 , A ,tanf3 
OgtPi) = m Si (R)-Mm^(pi)-97) 
In this case we have to add two more terms to Eq.(44), namely : 



(46) 



x 2 ^x 2 +(^^r+(^-^r (47) 



and minimize this expression by varying also 9g and 6^ . As one can see the errors on 9 g 
and Or do not contribute to the errors of the other parameters in the first approximation. 



As an example, in the Figs. [40] 4- 43 we have represented the dependence of the observ- 
ables on mo, mi/2, Aq and tan/3, respectively, in point 5. Each histogram is fitted with a 
polinomial in a domain ±50 GeV - for mo, ^1/2 an d Aq and ±0.5 - for tan/3, around the 
nominal values of point 5. The first derivatives of these functions are used to obtain A^ 
ofEq.( _ 



The final results are compiled in the form of relative errors in Table 11, taking Aq fixed 
at its nominal value of each point. The obtained precision on the SUGRA parameters are 
in general higher than in the case of conserved R parity ||. The reason is that here we 
were able to determine more observables and usually in a more direct way through the 
reconstruction of the LSP. 
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Table 6: The slopes, d P a , of the observables at the nominal value of point 5 : mo = 100 GeV. 
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Figure 41: The dependence of m-£«, m i R , m q R , m q L i m },a and M^ggp on smoothed by a 

polynomial fit. All other parameters are fixed at those of point 5. 



Table 7: The slopes, 
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Figure 42: The dependence of m^o, m i R > m qm m qL^ m h° anc ^ ^ossf on smoothed by a 
polynomial fit. All other parameters are fixed at those of point 5. 



Table 8: The slopes, d -f , of the observables at the nominal value of point 5 : Aq = 300 GeV. 
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Table 9: The slopes, 9 , ° a , of the observables at the nominal value of point 5 : tan/3 = 2.1 . 
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Table 10: The measured values of observables in SUGRA points 1, 3 and 5 for each ft coupling 
considered. All the values are in GeV. 
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Table 11: The relative errors on m , m 1 / 2 and tan/3 in SUGRA points 1, 3 and 5 for each ft 
coupling considered. 
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6 Conclusions 



We have studied the feasibility to detect a SUSY signal by ATLAS in the framework of 
the SUGRA model and to determine its parameters in the case when R parity is broken in 
conjunction with lepton number violation: A 7^ 0. For this purpose we have chosen three 
representativ points in the SUGRA parameter space and two different type of couplings, 
both having a value 10~ 3 , small enough to concentrate the effect in the LSP decay but 
large enough not to see displaced vertex in this decay. 
Our conclusions are the following: 

1. The SUGRA signal is visible in a very large domain of the parameter space, even 
beyond tuq ~ m^j ~ 1 TeV. 

2. The energy scale of the SUGRA signal can be determined by inclusive measurements 
like the effective mass (m e ff) or the normalised transverse momentum per lepton (p l t ' norm ). 

3. In the case of couplings with absence of a r among the decay products of the Xi 
(e.g. A122) one can reconstruct the SUSY particles and this reconstruction can be used for 
a precision determination of the model parameters. The achieved precision turns out to be 
better than it was the case with conserved R parity. This is because one can reconstruct 
the LSP from its decay products. At the low energy point where the chargino or second 
lightest neutralino produces additional leptons this determination is slightly handicapped 
by the combinatorial background and the most complex structure. 

4. In the case of LSP decay with a r particle in the final state the full reconstruction 
of the LSP, i.e. the determination of its four momentum, is not always possible, however, 
one can still estimate its mass (except at point 3 - A123). It allows ones to determine the 
parameters of the SUGRA model in spite of the large combinatorial background due to 
the leptonic decay of the LSP. This determination in most of the cases is better or at least 
comparable in precision with that when R parity is conserved. 
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